Explosive percolations on the Bethe Lattice 
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Based on the self-consistent equations of the order parameter P^, and the mean cluster size S, 
we develop a novel self-consistent simulation (SCS) method for arbitrary percolation on the Bethe 
lattice (infinite homogeneous Cayley tree). By applying SCS to the well-known percolation models, 
random bond percolation and bootstrap percolation, we obtain prototype functions for continuous 
and discontinuous phase transitions. By comparing the key functions obtained from SCSs for the 
Achlioptas processes (APs) with a product rule and a sum rule to the prototype functions, we show 
that the percolation transition of AP models on the Bethe lattice is continuous regardless of details 
of growth rules. 



PACS numbers: 64.60.ah, 64.60. De, 05.70.Fh, 64.60. Bd 



Since Achlioptas. et al. [l[ suggested an explosive per- 
colation model, there have been intensive studies on the 
explosive percolations (EP) @-[I2|- Achlioptas process 
(AP) was originally argued to show the discontinuous 
phase transition on the complete graph (CG) by sup- 
pressing growth of large clusters [l[ . Subsequent studies 
on variants of EP models on networks and lattices also 
argued to show the discontinuous transition [Ij-Q . In con- 
trast Riordan and Warnke [l3[ analytically showed that 
the phase transition in AP model [l| on CG is continuous 
by use of the arbitrary connectivity of CG. Furthermore 
several studies also showed that the transition of variants 
of EP models on CG is continuous [IcHl2l |. However, 
EP on CG was still reported to undergo discontinuous 
transition depending on the detail of cluster growth rule 
[HI, [l5j]. Therefore the transition nature of EP on CG is 
not still clear. Since the dimensionality of CG is infinite, 
physics on CG must satisfy the mean-field theory. In this 
sense the mean-field theory of explosive percolation is not 
still clearly understood. 

The Bethe lattice (infinite homogeneous Cayley tree) 
is physically a very important substrate or medium on 
which the mean-field theories for various physical models 
become exact (l6| . The analytic treatments of magnetic 
models \17\ . percolation [l6], HH and localization [l6| on 
the Bethe lattice give important physical insights to the 
subsequent developments of the corresponding research 
fields. One of theoretical merits of the Bethe lattice is 
that one can setup the exact self-consistent equations 
(SCEs). In this letter, by use of the exact SCEs we de- 
velop a novel self-consistent simulation (SCS) method for 
arbitrary percolation process on the Bethe lattice. From 
SCS method, we precisely calculate the order parameter 
Poo and the average size S of finite clusters of the AP 
models with a product rule or a sum rule. The obtained 
Poo and S can clarify the transition nature of AP models 



in the infinite dimension exactly. Furthermore unlike AP 
on CG, the bond connections on the Bethe lattice are 
purely local. Since there have been some papers that EP 
models on lattices with local bond connections show the 
discontinuous transition 0,0,11, it is physically impor- 
tant to study EP models on the Bethe lattice or in the 
mean-field level with local connections. As we shall see, 
the transition of AP on the Bethe lattice is continuous 
regardless of the sum rule or the product rule. This result 
physically means that AP models with local connections 
in the mean-field level shows the continuous transition 
irrelevant to the details of growth rules. In this sense 
the study on the Bethe lattice should give new physical 
insights to the mean-field theory of AP. 
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FIG. 1: Schematic diagram for a Bethe lattice with z — 3. 
The center part of the Bethe lattice is a 2-generation Cayley 
tree with 6 edge sites denoted by "X". Each edge site is 
connected to an infinite branch. Each infinite branch has an 
infinite cluster (IC) with the probability 1 — R(— A). An 
infinite branch has a finite cluster (FC) of the average size Sb 
with the probability R. 

From the structure of the Bethe lattice, SCEs for 
macroscopically measurable physical quantity can be eas- 
ily setup [l6| . Let's first setup SCEs for arbitrary perco- 
lation on the Bethe lattice. To setup SCE on the Bethe 
lattice in Fig. [T] one starts with a center site (or origin 
O) having z bonds. Then let's think about a part of the 
Bethe lattice with n generations from O, which have to- 
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tal N = 1 + z(k n - l)/(k - 1) sites with k = z - 1. To 
make complete Bethe lattice, one should add an infinite 
branch to each of zfc" -1 edge sites. To find the proba- 
bility that O belongs to the infinite cluster, we need to 
know the occupation probability p of a bond (or a site) 
and the probability R with which an edge site does not 
connected to the infinity through the infinite branch con- 
nected to it or the branch connected to the edge site has 
only finite clusters. Then the order parameter of an 
arbitrary percolation, which is defined by the probability 
of O to belong to the infinite cluster, is a function of n, 
p and R as Poo(p, R, n). Then SCE for is 

Poo = P 00 (p,R,n) = P 00 (p,iJ,n') (1) 

for any combination of {n, n'}. For the random or 
normal site percolation Eq. (TTJ) with the combination 
{n = 1, n' = 2} gives R = 1 - p + pR k [IE El, which an- 
alytically reproduces the mean-field properties of percola- 
tion transition. If the probability that a cluster including 
O with s sites and t edge sites within the n-generation 
tree is P(p, s, t, n, R, Sb), then 

Poo(p,R,n) = l-^i2 t ^P(p,s,t ) n,i?,5 6 ), (2) 

t s 

where Sb is the average size of finite clusters connected 
to an edge site of the n-generation tree in Fig. [TJ SCE 
for the average size S of the finite clusters including O 
can also be written as 

S = S(p,R,Sb,n) = S(p,R,Sb,n'). (3) 

S(j>, R, Sb, n), where 

„, D c , E a ^P(P,s,t,n,R,S b )[8 + tS b ]Ii* 
S(p,R,Sb,n) = — — . (4) 

If one cannot calculate P(p, s, t, n, R, Sb) analytically 
or if one need to know R and Sb apriori to occupy a 
new bond (or site), one should estimate P(p, s, t, n, R, Sb) 
indirectly to solve SCEs (fTJ) and ([3]). One of the in- 
direct methods is a simulation method. In this let- 
ter we develop a simulation method to solve SCEs, 
which wc call the self-consistent simulation (SCS). In 
SCS, P(p, s,t,n, R, Sb) is estimated by the relation 
P(p,s,t,n,R,Sb) = N(p,s,t,n,R,S b )/N c i uster , where 
N(p, s, t, n, R, Sb) is the number of clusters including O 
with s sites and t edge sites within the n-generation tree 
occurred in the simulation runs. Of course N c i uste r is 
the total number of clusters which includes O within 
the n-generation tree occurred in the same simulation 
runs. In the simulation both P(p, s,t,n, R, S b ) and 
P(p, s, t, n', R, Sb) are estimated simultaneously using the 
Bethe lattice with the center n-generation tree if n > n'. 
Since we don't know R(p) and S b (p) apriori, the itera- 
tion processes are needed in SCS. From initially guessed 
values for R{p) and S b (p), the final or saturated values of 
R(p) and S b (p) are obtained by the iteration of unit sim- 
ulation process. The unit simulation process consists of 



the following two steps I) and II). I) By use of the sim- 
ulation runs based on the given values R(p) and S b (p) 
P{p,s,t,n,R,Sb) and P{p,s,t,n' ,R,Sb) are estimated. 
II) From the estimated P{p, s, t, n, R, Sb), new R(p) and 
Sb{p) are calculated by utilizing SCEs (TTJ and ([3]). In the 
unit simulation process to get the new R(p) and Sb(p), 
the quantities like P(p, s, t, n, R, Sb) are estimated by av- 
eraging over at least 10 6 simulation runs. Such unit pro- 
cess is repeated until R{p) and S b (j>) reach the saturation 
values. Using the saturated values of R(p) and Sb{p), Poo 
and S are estimated from Eqs. © and (U]). 
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FIG. 2: (color online) SCS results of f p (A) and Poo for 
random bond percolation (RBP), bootstrap site percolation 
(BSP) and AP on the Bethe lattice with z = 4. (a) f p (A) 
of RBP for p = 0.31333(< p c ), p = 0.33333(~ p c = 1/3), 
p = 0.34333(> p c ) and p = 0.35333(> p c ). (b) Plot of P x for 
RBP against p. The line denote the analytic result, (c) f p (A) 
of BSP for p = 0.87888(< p c ), p = 0.88888(~ p c = 8/9), 
p = 0.89388(> p c ) and p = 0.89888(> p c ). (d) Plot of 
for BSP against p. The line denote the analytic result. In- 
set shows the evolution of results of iteration processes. The 
arrow — » denotes the direction from the earlier iteration pro- 
cess to the later iteration process, (e) f P (A) of AP with a 
product rule (PR) for p = 0.60575(< p c ), p = 0.61575(~ p c ), 
p = 0.61775(> p c ) and p = 0.62575(> p c ). (f) Plot of for 
AP with PR against p. Inset shows the evolution of results of 
iteration processes. The arrow — > means the same thing as in 
(d). 

The phase transition of random or normal percolation 
(RBP) on the Bethe lattice is well-known to be contin- 
uous We have applied our SCS to random bond 
percolation on the Bethe lattice with z = 4. The results 
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for z — 4, n — 13, n' = 5 arc displayed in Figs. [2ja) and 
(b). In Fig. [H[a) we display the simulation results for 



/pCAJsPoofel-^nO-Poofel-^n), (5) 

where A = 1 — P for various p. From Eq. ([1]), A* 
(or P*) which satisfies f p (A*) = is the real physical 
value for a given occupation probability p. For p < p c 
there occurs only trivial solution A* = as shown in 
Fig. [2ja). Increasing p from p c , the nontrivial solu- 
tion A* > continuously increases from zero. This con- 
tinuous increase makes the order parameter Poo(p) in- 
crease continuously as in Fig. [5{b). The simulation re- 
sult for Poo(p) exactly coincides with the analytic result 
Poo = 1 - [1/2 - ^(4 - 3p)/4p] 4 for p > p c (= l/k = 1/3) 
as show in Fig. [2jb). Therefore f p {A) which behaves like 
in Fig. [Ha) is a prototype function for the continuous 
transition. 

Bootstrap site percolation (BSP) on the Bethe lattice 
[l~9j is analytically known to show the discontinuous tran- 
sition. In BSP the order parameter Poo is the probability 
for an occupied site to be a site of the infinite m-cluster. 
Here the m-cluster means the cluster in which every oc- 
cupied site has at least m occupied nearest neighbors. 
For m > 3 the phase transition of BSP on the Bethe 
lattice is discontinuous (Isj . By using SCS we have ob- 
tained f p (A) and Px, for BSP with m = 3 on the Bethe 
lattice with z = 4. The results of SCS for BSP with 
z = 4, m = 3, n — 13, n' = 1 arc depicted in Figs. [2jc) 
and (d). As shown in Fig. EJ^c), f P (A) — for p < p c 
has only the trivial solution A* = as in RBP. In con- 
trast to RBP the nontrivial solution of f p (A) = for 
p > p c comes from the peculiar behavior of f p (A), which 
reminds us the thermodynamic instability in the thermal 
mean-field first order transition, which makes the sudden 
jump of the A* from zero at p = p c = 8/9. The jump 
of A* causes the discontinuous increase of Px, as shown 
in Fig. 2(d), which is exactly the same as the analytic 
result, Poo = p(l - P) 4 + 4pR(l - P) 3 with R = 1 for 
p < p c and R = [1 - 3^/1-^/^/4 for p > p c [ig[. 
Therefore f p (A) like in Fig. 2(c) is a prototype function 
for the discontinuous transition. Furthermore, as shown 
in the inset of Fig. 2(d), the iteration processes in SCS 
for the BSP with the initial value R = drives Poo(p) 
from the continuous increase to the final discontinuous 
jump at p = Pc, which is also a typical behavior of the 
discontinuous transition on the Bethe lattice. 

We now focus the AP model. To occupy a bond in 
the AP model [l[, two bonds a and /? are randomly 
chosen. Let S a i and ^2(5^1 and S/32) be the sizes of 
the two clusters which would be connected by occupying 
bond a (/3). Under a product rule (PR), the bond a is 
chosen for the growth of clusters and f3 is discarded if 
llj=i Saj < Spj- Otherwise, the bond (3 is cho- 

sen. AP model with a sum rule (SR) is the same as 
that with PR except for the change of the condition into 
E S a j < Sf)j- Note that arbitrary edge site is 



connected to an infinite cluster with the probability A 
depicted as in Fig. Q] Therefore S a (p)j is calculated as 

Sa(/3)j = s a(P)j + OO X I a (/3)j + (*as(/3)j - I a (P)j)Sb , (6) 

where s Q (^)j is the number of sites within the n- 
generation tree in the cluster a(/3)j, t a (mj is the num- 
ber of edge sites and I a (/3)j is the number of edge sites 
which are connected to the infinite cluster. Of course 
I a {fi)j depends on R or A. Therefore P(p, s, t, n, R, Sb){= 
N(p, s, t, n, R, Sb) / N duster) depends apriori on R and Sb 
and thus the iteration is essential in SCS for AP model. 
For SCS of AP model, it should be careful to choose 
n'(< n) for a given n, because too small n' cannot con- 
vey enough information for AP and n' close to n can 
never gives physically plausible solutions for SCEs ([T]) 
and ([3]). From the simulations with various {n,n'} it is 
confirmed that the suitable choice of n' should be in the 
interval n/3 < n' < n/2. f p (A) and P^ for AP model 
with PR obtained from SCS with z = 4, n = 14, n' = 5 
are displayed In Figs. [He) and (f). The results for AP 
model with SR are nearly the same as those in Figs. (2Je) 
and (f) except that p c for SR is slightly smaller than p c 
for PR. As can be seen from Fig. [2je), f p (A) = has 
only trivial solution A* = for p < p c as RBP and BSP. 
Increasing p fromp c , the nontrivial solution A*(> 0) con- 
tinuously increase from zero. A* for AP increases very 
rapidly compared to A* for RBP as p increases. Except 
for this rapid increase, f p (A) for AP behaves in the same 
way as that for RBP in Fig. [H[a), which is the prototype 
function for the continuous phase transition. We cannot 
find any symptom that f p (A) for AP behaves like f P (A) 
for BSP in Fig. 2(c). The continuous increase of A* 
makes the order parameter Pca(p) for AP increase con- 
tinuously as in Fig. HJf). Moreover, as can be seen from 
the inset of Fig. [5Jf) the iteration processes in SCS for 
AP with PR drives Poo(p) from the discontinuous jump 
to the final continuous increase at p = p c , contrary to 
those in the inset of Fig. [2jd). 
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FIG. 3: Poo 1 against Ap(= 1/(N - 1)) for {n,n'} = 
{7, 4}, {9, 4}, {11, 5}, {13, 5}, {14, 5}. (a) AP model with PR 
and (b) the model with SR (b). 

In SCS on the Bethe lattice with the n-generation cen- 
ter tree, we cannot obtain P^ against all the continuous 
p G [0,1]. Instead in SCS p is increased discretely by 
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Ap = 1/(N — 1), where N is the total number of sites 
in the n-generation Cayley tree. Therefore in SCS the 
lowest nonzero occurs at p — p' c very close to the 
true p c with < (p' —p c ) < Ap. In Fig. [3]we display the 
dependence of Px, on Ap for SCSs of AP models with 
{n,n'} = {7, 4}, {9, 4}, {11, 5}, {13, 5}, {14, 5}. As can 
be seen from Fig. O P^ for both AP models with PR 
and SR decreases monotonically to zero as Ap decreases 
to zero. This result for also supports the fact that 
the phase transition nature of AP on the Bethe lattice is 
continuous. 
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FIG. 4: S on the Bethe lattice with 2 = 4. (a) Plot of S for 
RBP against p. The data are obtained from SCS with n = 13, 
n = 5. (b) Plot of S for AP with PR against p. The data 
are obtained from SCS with n = 14, n' = 5. 

From SCS based on Eq. ([3]), the average size S of the 
finite clusters is obtained for RBP and AP model with 
PR as in Fig. [4] S for AP model with SR shows nearly 
the same behavior as Fig. Hl[b). It is confirmed that S for 



RBP from SCS is nearly identical to the analytic result 
5 = 1+ 4pR/(l - p - 2pR), where R = 1 for p < p c and 
R=l + ^/(4 - 3p)/4p 3 - 3/2p for p > p c . Even though 
S for AP model diverges more rapidly than S for RBP, 
S for both RBP and AP diverges as S(p) ~ \p — p c |~ 7 
with the susceptibility exponent 7 = 1.00(1). The result 
for S of AP models also supports that the transition in 
AP models is continuous on the Bethe lattice. 

Since the transition in the AP models is continuous, 
the order parameter exponent ft is calculated by fitting 
the relation P^ ~ (p — p c )^ to the data for p > p c but 
very close to p c in Fig. HJf). Since f p (A) for p ~ p c 
in Fig. [5{e) increases very slowly as A increases from 
zero, Poo increases very rapidly and the exponent ft is 
expected to be very small. From the best fit the obtained 
exponent is ft — 0.05(5). The data even fits very well to 
Poo ~ I ln(p — p c )\~ x with x = 3.4(1). This result means 
that the exponent ft on the Bethe lattice is very small 
and nearly identical to those obtained on the complete 
graph 0,0, [3. 

In summary, we show that AP models on the Bethe 
lattice have a continuous transition from the SCSs de- 
veloped to cover arbitrary percolation processes on the 
Bethe lattice. For this f P {A) for AP is first shown to be 
physically identical to that for RBP. We also shown that 
Poo 1 for AP model decrease to zero as Ap = l/(N — 1) 
goes to zero. The divergent behavior of S for AP mod- 
els is also shown to be the same as that for RBP with 
7=1. The exponent ft is also shown to be very small or 
ft ~ 0.05. 
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